Quantum Characterization of a Werner-like Mixture 
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We introduce a Werner-like mixture [R. F. Werner, Phys. Rev. A 40, 4277 (1989)] by considering 
two correlated but different degrees of freedom, one with discrete variables and the other with 
continuous variables. We evaluate the mixedness of this state, and its degree of entanglement 
establishing its usefulness for quantum information processing like quantum teleportation. Then, 
we provide its tomographic characterization. Finally, we show how such a mixture can be generated 
and measured in a trapped system like one electron in a Penning trap. 
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I. INTRODUCTION 



It is nowadays well known that the nonlocal properties 
of Quantum Mechanics enable striking processes in 
quantum information In these processes it is promi- 
nent the role of maximally entangled states |^J^. How- 
ever, very often, the decoherence effects due to the envi- 
ronment transform the pure entangled state into a sta- 
tistical mixture and degrade quantum entanglement in 
the real world Although purification schemes may 
be applied to noisy channels ||7|, there exist some mix- 
ture states which maintain interesting properties. An 
illuminating example is provided by the Werner mixture 
1^, which is not a mixture of product states, nonetheless 
not violating any Bell's inequality p], but still useful for 
quantum information processing [pyiO|]. Such states be- 
long to systems with two discrete degrees of freedom like 
two spin-i. However, information processing may some- 
times involve hybrid systems where one degree of freedom 
has discrete variables and the other continuous variables. 
It may happen, for instance, in trapped ions ]ri| , p^ , or 
in trapped electrons |13| or in cavity quantum electro- 
dynamics ]r^ , p^ . Thus, it will be the aim of this paper 
to consider a mixture, which resembles the Werner one, 
but one of the two subsystems is described by continuous 
variables. 

On the other hand, states and processes used in quan- 
tum information typically need of a well characterization 
1^ . This can be accomplished by using tomographic tech- 
niques . Concerning the quantum state measurement, 
after the seminal work by Vogel and Risken , a lot of 
progress has been obtained and further techniques and 
algorithms were developed p^ . We would just mention 
the possibility of state reconstruction, for a composite 
system of discrete and continuous variables, by simply 
measuring the set of rotated spin projections and dis- 
placed number operators, |l^JT^. Then, we shall provide 
the tomographic characterization of a Werner-like mix- 
ture by generalizing that method. 

The outline of the paper is the following: In Section II 
we discuss the Werner mixture and we extend the concept 
by considering one of the two subsystems with continu- 



ous variable. Then wc characterize such a state in terms 
of mixedness and entanglement. Section III is devoted 
to the tomographic method employed for such a state 
reconstruction. In Section IV we present the results of 
numerical simulations. Finally, in Section V we discuss 
a possible implementation and Section VI is devoted to 
conclusions. 



II. WERNER-LIKE MIXTURE 

In his pioneering paper. Bell proved that a local real- 
istic interpretation of Quantum Mechanics is impossible 
1^, and for the case of pure states it is known that, when 
measurements are performed on two quantum systems 
separated in space, their results are correlated in a man- 
ner which, in general, cannot be explained by a local 
hidden variables model ||2^. Since the only pure states 
satisfying the Bell inequality are pure product states, 
one might naively think that the only mixed states that 
do not violate Bell's inequality are mixtures of product 
states. However, Werner showed that this conjecture 
is false for the so-called Werner states 
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where li [i — 1,2) stands for the identity operator of 
a single qubit (| | + | t>«(T I) and 1^--) = 

^ (I T)2 - I T)i| 1)2) is the spin singlet state. 

A more general Werner mixture can be obtained by 
considering one of the two subsystems, say 2, as described 
by continuous variables. A way to encode qubit in con- 
tinuous variable systems could be the use of even and odd 
cat states which are orthogonal |^l| , thus resulting in the 
same situation of Eq.([l]). Instead, the choice we are go- 
ing to make is more general and gives the possibility to 
explore a variety of situations. 

That is, we now replace the states | 1)2, | i)2 of the 
second qubit with \a)2 and | — a)2, where the latter are 
coherent states of amplitude a and —a respectively (we 
shall consider a S R throughout the paper for the sake 
of simplicity). Therefore, a Werner- like mixture would 
be 
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Since (a| — a) ~ exp(— 2q;^) 7^ 0, the above state doesn't 
describe a real two qubit system, but rather a two qubit 
system with nonorthogonal states |2^]. Of course, for 
a ^ 1 Eq.(||) behaves hke the state (|^), but we want 
to study its characteristics for a generic value of a. To 
this end, we map the above state in the two spin-^ 
Hilbert spaces by introducing for the subsystem 2 a vec- 
tor 1-0)2 = k\ 1)2 +vl — f?\ t)2 and considering the non 
orthogonal states | 1)2 , |V-')2 instead of | — a)2 , \c()2, with 
("01 I) — 1^ = exp(— 2a^). Then, Eq.(^ can be rewritten 
as 

p-^{u>iiai®u)22ai 

+1 T)ii(T l®U)22(i I 

+ 1 |^|V')22(V'l 
+ 1 T)ll(T l®IV')22(V'|} 



+ ^[| i)l|V^)2-|T)lU)2] 
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The degree of mixedness of the state (|) can be evaluated 
by using the Von Neumann entropy [ 231 

S{p) = -Tr (plog2 p) = - ^ A, log2 A, , (4) 



where Xi are the eigenvalues of the matrix representation 
of p. We have calculated such eigenvalues in the basis 
B={\ i)2, I T)2, I T)i| i)2, I T)i| 1)2} and we 
have plotted the entropy S in Figj^. 
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We can see that also for a = the state is a mixture, 
then, it becomes more and more mixed by increasing the 
value of a, but never reaching the maximum (a com- 
pletely mixed density operator in a d-dimensional space 
has entropy log, d). It is also worth noting that for a (k) 
arbitrary, Eq. (J^) is not a mixture of Bell's states as it is 
for the Werner mixture (|]), i.e., for a ^ 1. 

Let us now consider the measure of entanglement. For 
a two spin-^ system the state separability is related to 
the partial transposition operation p.p^ . The matrix el- 
ements of partial transposition p'^^ of a state p are given 



by p 
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where 
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with TO, n, /i, v =1, |. A density matrix p for a two spin-2 
system is inseparable if and only if its partial transpose, 



measure of entanglement can be defined as [Eq 



p^, has any negative eigenvalue [jipj]. Then, a suitable 



(5) 



E{p) = -2J2K 



where A~ is a negative eigenvalue of . It is worth not- 
ing that this measure satisfies the necessary conditions 
required for every measure of entanglement Then, 
we have plotted in Figj^ the quantity E calculated by 
exploiting again the matrix representation of Eq.(H) in 
the basis B. 
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FIG. 1. Entropy S as function of a = -y/ — (1/2) In k. 
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FIG . 2. Degre e of entanglement E as function of 
a = A/-(l/2)ln«:. 

We can see that the Werner-like mixture is factorisable 
for a = 0, then, by increasing a, its degree of entangle- 
ment increases saturating at the value 1/4 characteris- 
tic of the true Werner mixture (^. This latter value of 
entanglement is known to be sufficient to improve the 
teleportation performances over classical limit once the 
sender (Alice) and the receiver (Bob) initially share the 
state (^ Q. Thus, it is straightforward to ask to what 
extent a Werner-like mixture can be used for the same 
goal. To establish a threshold value for a we are going to 
consider the teleportation fidelity. To this end, we first 
write the state p through the Hilbcrt-Schmidt decompo- 
sition 
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where {crn}n=i,2,3 are the standard Pauh matrices, r, s 
are vectors in and r-cr — Y^^i=i f'i'^i- Furthermore, the 
coefRcients tn^m — Tr[p ®a'^] form the real matrix T 
describing the correlations between the two qubits. Thus, 
the teleportation capabilities will depend on the specific 
form of T. In particular it is shown in Ref. that the 
teleportation fidelity amounts to 
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Then, in Fig.^ we have plotted the quantity F versus a 
compared with the classical fidehty 2/3. 
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FIG. 3. Teleportation fidelity as 

function of a = y^— (1/2) In k. The dashed line represents 
the classical bound of 2/3. 

We immediately recognize the presence of a threshold 
value (a = 0.2476) below which the Werner-like mixture 
becomes useless for quantum teleportation. 



III. STATE MEASUREMENT 

We now discuss the possibility of a complete character- 
ization of the Werner-like mixture through tomographic 
techniques. In particular, we generalize the method pre- 
sented in Rcfs. Jl8| , p^ to non pure states. 

Accordingly to the state reconstruction principle devel- 
oped in Ref. we choose an observable, here ® a^a, 
then we apply suitable unitary transformations to get 
a set of observables giving the whole state information 
upon measurements. In our case the transformations 
would be 

U{d, If) — exp [—i9 ((Ti cos ip + 02 sin ip)] , (8) 
D{(3) = exp [(3a^ - /3*a] , (9) 



which lead to rotated (by angles and tp) spin projec- 
tion in the subsystem 1 |p8| , and to displaced number 
state (by a complex amount (3) in the subsystem 2 [ p9[ . 
Then, it is possible to consider the following measurable 
marginal distributions 

w{s,n-e,^,P) ^Tv{pD(P)U{e,^) 

\s)ii{s\®\n)22{n\U\e,p)D\P)) , 

(10) 

having as variables the eigenvalues s, n of (73 and a^a 
as variables and parametrically depending oti 9^ p and 
f3. Thus, measuring the state p would mean the possi- 
bility to express p as a functional operator of w, i.e. to 
invert expression (^0|). This also means the possibility to 
sample the density matrix elements (in some basis) from 
the quantity measured by spanning the whole space of 
parameters. 

In reality, we shall see that it is not necessary to con- 
sider all possible values of parameters. As matter of fact, 
we write the density operator (H) as 
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(11) 



where each operator p^'^ ^ p^^, pJ-^, p-^^ , can in turn be 
represented in the Fock basis of the subsystem 2. 

Now, we set 6 = p — Q and we suppose to retain 
only the measurement results s —]] then, expanding the 
density operator p^^ in the Fock basis, and defining Nc 
as an appropriate estimate of the maximum number of 
excitations (cut-off), we have 



u;(T,n;0,0,/3) 



E 

fc,m— 



[n,(3\k){k\p^''\m){m\n,(3) . (12) 



The projection of the displaced number state |n, (3) onto 
the Fock state \m) can be obtained generalizing the result 
derived in Ref. Q. 

Let us now consider, for a given value of |/?|, w as a 
function oi (f) — arg[/3] |Q and calculate the coefficients 
of the Fourier expansion 



w(t, n; 0, 0, r) 



27r 



d0u;(t,n;O,O,0)e''''^ , (13) 



for r = 0, 1, 2, . . .. Combining Eqs. (|T2|) and (|l|), we get 

Na-r 

z/;(T,n;0,0,r)= ^ GW„(|/3|)(m + r|pTT |^) , (14) 

m=0 

where the explicit expression of the matrices G is given 
in Ref. |l|]. 

We may now notice that if the distribution w{'\ 
, n; 0, 0, 0) is measured for n S [0, N] with N > then 
Eq. ([1^) represents, for each value of r, a system of A^-|- 1 
linear equations between the + 1 measured quantities 
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and the Nc + 1 — r unknown density matrix elements. 
Therefore, in order to obtain the latter, we only need to 
invert the system 



N 

(m + r|p"|m) = ^MW„(|/3|)«;(T,7i;0,0,r), (15) 

n=0 



where the matrices M are given by M = (G'^G) ^G^ . 

The procedure can be repeated with w{l, n; 0, 0, /?) in 
order to get the matrix elements of p^^. Then, chang- 
ing the parameters so that 9 — tt/4 and (p = — 7r/2 we 
can get the real part of the matrix elements of p'^K In- 
stead, with 9 = IT I A and Lp = Q we can get the imaginary 
part of the matrix elements of p^-'-, thus concluding the 
reconstruction procedure of the state (11). 



IV. NUMERICAL RESULTS 



As an example of the proposed method, we show in 
Figs. ^, 1^ the results of numerical Monte-Carlo simu- 
lations of the reconstruction of the state (||) once written 
in the form (^Tj). In this simulation we have used the 
value a = 0.7 which make the state different from a true 
Werner mixture, but still having nonclassical features as 
discussed in Section II. 

In order to account for experimental conditions, we 
have also considered the effects of a non-unit quantum 
efficiency rj in the counting of the number of excitations. 
When 77 < 1 , the actually measured distribution is related 
to the ideal distribution by a binomial convolution [^2| . 

Statistical errors are accounted for as well by consid- 
ering an estimation of the marginal distributions given 
by w^^^{n](f)) = N,j,{n)/M^^ where N^{n) is the number 
of events with n counts at phase <j), while A/^ is the total 
number of events at the same phase. Then, following the 
arguments given in Refs. pl| , p3| , the quantities Af^{n) 
can approximately be regarded as independent Poisso- 
nian random variables, whose means and variances are 
given by w{n; (j))/^^ « w'^^^{n;(j))J\f^. The variance of 
w'^^^ {n; (/)) may then be approximated by ■w''^^{n;(j))/J\f^, 
so that the variances of the real and imaginary parts of 
the density matrix can be easily estimated using Eq. ( p^ ) 
(and similarly for the other density operators). This 
means that the errors can be estimated in real time as the 
experiment runs, simultaneously with the reconstruction 
of the density matrix elements. 

Other error sources leading to discrepancies between 
true and reconstructed density matrices can be identified 
in the choice of \(3\, Nc and A^^ (number of phases). How- 
ever, as it can be seen from Figs. ^, H, |^ the reconstructed 
density matrices turn out to be quite faithful. 







n \r n 
FIG. 4. Density matrix elements of p^^ for the state 
with Q = 0.7. In the computer simulation — 96 phases 
are considered and A/0 = 10* measurement events at each 
phase are assumed to be recorded using Nc — 31. The other 
parameters are |/3| = 0.6 and rj = 0.9. (a) Exact density 
matrix; (b) reconstructed density matrix; (c) statistical er- 
rors; (d) absolute difference between reconstructed and exact 
density matrix elements. 







FIG. 5. Density matrix elements of p^^ for the state 
with a = 0.7. In the computer simulation N^ — 96 phases 
are considered and Afev ~ 10* measurement events at each 
phase are assumed to be recorded using Nc — 31. The other 
parameters are \(3\ = 0.6 and rj = 0.9. (a) Exact density 
matrix; (b) reconstructed density matrix; (c) statistical er- 
rors; (d) absolute difference between reconstructed and exact 
density matrix elements. 
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FIG. 6. Density matrix elements of p^^ for the state (^) 
with a = 0.7. In the computer simulation A^^ = 96 phases 
are considered and JVev = 10* measurement events at each 
phase are assumed to be recorded using Nc — 31. The other 
parameters are |/3| = 0.6 and rj = 0.9. (a) Exact density 
matrix; (b) reconstructed density matrix; (c) statistical er- 
rors; (d) absolute difference between reconstructed and exact 
density matrix elements. 



In addition to the present case we have performed other 
simulations with different values of a and several values of 
the parameters, which confirm that the present method 
is quite stable and accurate. 



The phase-space description corresponding to 
given by the Wigner-function matrix |34] 



IS 



W = 



Tr2[pTT<5(^_a)] Tr2[pT^ ^(7 - a) 
Tr2[p" 6{-f - a)] Ti2[p^^ S{-f - a) 



(16) 



where ^(7 — a) is the Fourier transform of the displace- 
ment operator ||30|l . 

Then, the Wigner functions corresponding to the den- 
sity matrices of Figs.^, || and ^, are shown in Figs.^ and 
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FIG. 7. Top: Wigner function corresponding to the den- 
sity matrix of Fig.^ (true on the left and reconstructed on the 
right). Bottom: Wigner function corresponding to the den- 
sity matrix of Fig.^ (true on the left and reconstructed on the 
right). 
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FIG. 8. Wigner function corresponding to the density ma- 
trix of Fig.0 (true on the left and reconstructed on the right) . 



The reconstructed Wigner functions as well turn out to 
be quite faithful. We would like to emphasize their par- 
ticular shape. In W^^ there are two hills, one centered 
in Re 7 = —a and the other in Re 7 = a coming from 
the random part of the state (§), instead, the pseudo sin- 
glet part of the state contributes only to the bump in —a, 
thus producing the asymmetric effect. The opposite hap- 
pens for W^K The shape of W^i is due to the quantum 
interference given by the entanglement between the two 
degrees of freedom: in fact, in absence of entanglement 
p^^ would just be a replica of the diagonal parts p^T and 
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V. PHYSICAL REALIZATION 

We now briefly discuss a system where the Werner-like 
mixture could be synthesized and even measured. It is a 
single electron trapped in a Penning trap |^5| , where two 
different degrees of freedom of the same particle can be 
entangled and also measured. 

An electron in a Penning trap is confined by the combi- 
nation of a homogeneous magnetic field along the positive 
z axis and an electrostatic quadrupole potential in the xy 
plane . The spatial part of the electronic wave func- 
tion consists of three degrees of freedom, neglecting the 
slow magnetron motion (whose characteristic frequency 
lies in the kHz region), here we only consider the axial 
and cyclotron motions, which are two harmonic oscilla- 
tors radiating in the MHz and GHz regions, respectively. 
On the other hand, the spin dynamics results from the 
interaction between the magnetic moment of the electron 
and the static magnetic field, so that the free Hamilto- 
nian reads as 135(1 



^^frcc = 'hujzO^^az + tkiJcaluc + TujJsaz/2 . 



(17) 



where the indices z, c, and s refer to the axial, cyclotron 
and spin motions, respectively. 

Here, in addition to the usual trapping fields, we con- 
sider an external radiation field as a standing wave along 
the z direction and rotating, i.e. circularly polarized, in 
the xy plane with frequency fl To be more specific, 
we consider a standing wave within the cylindrical cavity 
configuration p7[ | with the (dimensionless) wave vector 
k. Then, the interaction Hamiltonian reads pq] 



Hint = ?ie 



cos(fcz + x) 



+ K L e*("*+'^) + cr+e-'("*+^)j sin(/cz -|- x) , (18) 

where a± — (cri ± icr2)/2, and z = az + a\. The phase x 
defines the position of the center of the axial motion with 
respect to the wave. Depending on its value the electron 
can be positioned in any place between a node (x = 0) 
and an antinode (x = i'l'/S). Instead, the phase tp is 
related to the initial direction of the electric (magnetic) 
field in the xy plane, or to the phase of any reference field. 
The two coupling constants e and C a-re proportional to 
the amplitude of the applied radiation field. Depending 
on VL and Xi the interaction Hamiltonian (^8|) gives rise 
to different contributions at leading order in the Taylor 
expansion of sin(/cz -I- x) and cos(A;2; -I- x)- 

We immediately recognize the possibility of imple- 
menting the transformations (H), (^ on the spin and 
cyclotron degrees of freedom by appropriately exploiting 
the Hamiltonian (p^). For instance, U can be realized 
by setting 17 = w^, X = Oj a-iid then adjusting (p and 
= C^t. Differently, D can be realized by setting Vl ~ u>c, 
X = and then adjusting /3 = ete~"^. 

These transformations easily allow to generate the dis- 
entangled components of the mixture (||), i.e., | t)i|±a)2 



and I ± a)2, starting from the typical initial state 

IT)i|0)2. 

Instead, for what concerns the generation of the entan- 
gled fraction of the mixture (H), we recall the procedure 
developed in Rcf. That is, we consider the possi- 

bility of introducing pulsed standing waves through the 
microwave inlet psf so that e, C become time dependent 
and e, <^ indicate the pulse area (the duration of the pulse 
is assumed to be much shorter than the characteristic ax- 
ial period, which is of the order of microseconds). Then, 
nonclassical cyclotron states can be entangled with the 
spin states through the following steps . 

• First, we consider x = 0, = and a pulsed 
standing wave lasting Ati — t\ — = t\\ 

• Second, we allow a free evolution for a time At2 = 
t2 - ti = V(2c^^); 

• Third, we consider the action of another pulsed 
standing wave with x — ^■?r/2, = Wc, for a time 

Finally, if we consider the initial axial state as a Gaus- 
sian state with momentum width much smaller than C 
(which is easily obtained in the case of the ground state 
of the axial oscillator) , we end up with an evolution oper- 
ator of the form D(a.a\)^ where a is related to e, C, and p. 
It is then immediate to see that the initial state | t)i|0)2 
may evolve with the aid of D(aa\) and a spin rotation 
into 
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72 



(U)i|a)2-|T)i|-«)2) 



(19) 



This state has been already discussed in Refs. and 
constitutes the pseudo singlet component of the mixture 
(I). 

Thus, at each run of the experiment the desired com- 
ponent of the Werner-like mixture can be synthesized, 
thus allowing the generation of the state (||) on average 
ensemble. 

For what concerns the measurement, the addition of a 
particular inhomogeneous magnetic field (known as the 
magnetic bottle field |Q) to that already present in the 
trap, allows to perform a simultaneous measurement of 
both the spin and the cyclotron excitations number. The 
useful interaction Hamiltonian for the measurement pro- 



bottle 



— hujh 



■ t , 9 ' 



(20) 



where the angular frequency ujb is directly related to the 
strength of the magnetic bottle field. 

Equation (20) describes the fact that the axial angular 
frequency is affected both by the number of cyclotron ex- 
citations flt Oc and by the eigenvalue of . The modified 
(shifted) axial frequency can be experimentally measured 
|p5[ after the application of the inhomogeneous magnetic 
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bottle field. One immediately sees that it assumes dif- 
ferent values for every pair of eigenvalues of ajcc and 0-3 
due to the fact that the electron g factor is slightly (but 
measurably [^) different from 2. 

However, prior such kind of measurement, one has to 




which can be real- 
as discussed above. 



deal with the transformations ( 
ized through the Hamiltonian (] 

Repeating this procedure many times allows us to re- 
cover the desired marginal distributions, hence to sample 
the density matrix elements. 



VI. CONCLUSIONS 

In conclusions, we have studied the properties of a 
Werner-like mixture and a reliable method to achieve its 
tomographic characterization. A useful system to inves- 
tigate such states has been individuated in the trapped 
electron. There are other candidate systems that offer 
the possibility to generate and manipulate the studied 
state. We mention for example trapped ions or 
atoms in cavity quantum electrodynamics p^ . More- 
over, in such systems the studied state might involve two 
particles, or quite generally the two subsystems could be 
spatially separated. 

The experimental studies on this state might yield new 
insight in the foundations of quantum mechanics and al- 
low further progress in the field of quantum information. 
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